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ABSTRACT 

A study of the production and the decay of scalar and vector leptoquarks at 
high-energy e + e~ colliders is presented. All tree-order contributions to single 
leptoquark production have been calculated and incorporated in the Monte 
Carlo event generator ERATO-LQ. 
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Nature of physical problem: Leptoquark states emerge in several theoretical frameworks, 
including the Pati-Salam model, the MSSM, and composite models of quarks and leptons. 
Leptoquark states coupled to the first generation leptons and quarks can be studied in 
e + e~ colliders. The dominant mode, for leptoquark masses close to the energy threshold, 
is the single leptoquark production. Leptoquarks can be studied by selecting events with 
letpons and jets and by investigating the invariant mass distribution of the lepton plus 
jet system. 

Method of solution: Based on the four-fermion event generator ERATO, we have included 
all possible leptoquark contributions to the processes: 

e~e + — > £~£ + qq 

and 

e~e + — > £~ueqq' . 

All tree-order Feynman graphs as well as all necessary phase space mappings have been 
calculated and incorporated in the present extension ERATO-LQ. Moreover interference 
terms with the tree-order SM contributions have been fully included. The decay of the 
leptoquarks has been treated exactly. An additional code providing the total cross section 
in the resolved photon approximation has also been included. 
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1 Introduction 



Search for new particles has been constantly the center of attention in contemporary 
high energy physics. The initially reported excess of high Q 2 events at HERA by the 
HI 0] and ZEUS @ collaborations [] stimulated a substantial research interest f|, [|, || [7] 
on the possible existence of new particles of leptoquark type. Such particles emerge in 
several contexts, including the Pati-Salam model || and the grand-unified supersymmetric 
theories ||. Leptoquark- type interactions, described by higher- dimensional contact terms, 
are naturally incorporated within the idea of compositeness of quarks and leptons (|TIJ . 

It is the aim of this paper to present a complete description of the relevant four fermion 
final states, associated with leptoquark production at high-energy e + e~ colliders. This 
automatically includes pair leptoquark production as well as single leptoquark production 
channels. As a first application we present results concerning single leptoquark production 
at LEP2. 



We start the presentation with the most general effective Lagrangian ||16|| , which is 
SU(3) x SU(2) x U(l) invariant, baryon- and lepton-number conserving as well as family 
diagonal, in order to avoid the most severe constraints from low-energy considerations JT7J . 
It can be written as 

£f=-2 = (9ol q c L ir 2 iL + 9m u R e R ) S + g 0R d c R e R S + g 1L q c L ir 2 Ti L Si 

+ {g lL dl^h + gi R &fe R ) Vi„ + gi L u R Yh% + cc (1) 

£f=o = \hi L u R t L + hi R q~LiT2e R ) Si + h± L d R £ L Si + h R UR^e R V 0lJl 

+ (hoL qvfh + h 0R d R ^e R ) V Qfl + h XL q L =r^t L V Xii + cc (2) 



where F = 3B + L denotes the fermion number. The gauge-interaction terms are given 
by 

C = D^<$> - M 2 $ f $ (3) 

for scalars and 



£ = --Gl u G^ + M 2 $^ (4) 
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for vectors, with 

Gfiv = D^ v - D v ^ (5) 
and the SU(2)l x U(1)y covariant derivative given by 

Dfj, = [d„ - ieQ^Ap - %eQ z Z il - ieQ w (W+I + + W~r)] (6) 

where J ± are the 577(2) generators in the representation of the corresponding leptoquark 
state. The electroweak charges are given by 

cos 0,„ sitiL 



L For a recent review on leptoquark searches in high Q 2 events at HERA see reference ji) 
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and 

Q W = ■ (8) 

V2 sin 9 W 

In order to account for a more general structure in the vector leptoquark case one can 
also add the so called 'anomalous couplings' terms as for instance the 'magnetic' dipole 
one given by 

C extra = -i £ g V K v ^V^^ v (9) 

V=-y,Z 

parametrized in terms of k 7 and Kz- 



2 Single leptoquark production 



Leptoquark contributions to e + e~ processes can be classified into real and virtual ones. 
To the lowest order, virtual contributions proceed via the reaction e + e~ — > qq [|], §, [7], 
whereas for real ones we have the single leptoquark production mode |18|, [19|, |2(| and for 



light enough masses, or high enough energies, the pair leptoquark production |21 






Figure 1: EQ20, EQ10 and EQ5 classes. The curly lines can be 7, Z or W, with the 
appropriate changes in the fermion species. The other half of the graphs are obtained, as 
usually, by the interchange e{j>2) <-> e(j>3)- 

Leptoquark production at high-energy e + e~ collisions results to a four-fermion final 
state. We have calculated all tree-order contributions, including all interference terms 
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with the Standard Model, to the following 'semileptonic' reactions 



e~e+ -> £-£ + qq (10) 

and 

e"e + — > £~i> e qq , (11) 

where £ = e, fi,r. A full description of the MC generator ERATO can be found in ref- 
erences [p2] , p3f| . A brief description of its present extension is given in the Appendix. 
Leptoquark contributions to v?i>iqq and qqq'q' final states have not been considered here, 
since, as we will see, they are relatively suppressed. 

For the first reaction, Eq. ([TP) , with £ = e (£ = fi, r), the Standard Model Feynman 
graphs are of the NC48 (NC24) type0, whereas the leptoquark contribution can be classified 
in two categories: the LQ22 (LQll), which is diagrammatically identical to the Standard 
Model CC22 (CCll) with the W boson replaced by a leptoquark, and a new category, 
called EQ20, which contributes only when £ = e. The generic graphs belonging to the 
latter class are shown in Fig.|I[ Moreover there are nine scalar and nine vector leptoquark 
states contributing to these processes, namely 



qL qR qR qL qL qL qR qR J qR, 

and 

Vt + ,V^ + ,V^_,V^ + V L ,V R ,V R ,V^ + and , 

with the obvious notation Xg c , where a = L, R is the helicity of the lepton, b = I is the 
isospin of the leptoquark and c = I3 stands for its third component. 

For the reaction Eq.flllD, Standard Model contributions belong to the well known 
CC20(CC10) category, whereas those of leptoquarks fall into three different classes. The 
first one is the LQ20(LQ10), in close analogy to the previous case. The other two, the EQ10 
and EQ6, Fig.|l], receive contributions only from the first generation leptoquarks. In the 
L020(L010) class only Sq,Si and Uq,Ui states contribute, whereas in the classes EQ10 

and EQ6 we have contributions from S±, V± and S±, Vf 1 respectively. 

2 2 

In addition to the calculation of all tree-order Feynman graphs, we have also employed 



all phase-space mappings |24], [Eq| , which are necessary to cover all the kinematical regions 
where leptoquarks have the most substantial contributions. Initial state radiation (ISR) 
has also been included in the structure function approach [p5|] . Our calculation provides 



therefore the necessary framework to investigate real leptoquark production at high-energy 
e + e _ colliders, including both pair- and single-production modes, for all leptoquark types, 
as given by Eq.(|I|)-Eq.(^]), and for all generations. Nevertheless we find it convenient, to 
focus our present analysis on the first generation leptoquarks. It should be mentioned, 
however, that all necessary contributions for the study of the 2nd and 3rd generation 
leptoquarks at LEP2 have been fully accounted for in the present calculation. 

2 For the nomenclature about the four-fermion production diagrams see reference p3|. 
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(d) 



(e) 



Figure 2: The subset of Feynman graphs giving the dominant contribution to single 
leptoquark production: (a,b,c) signal, (e,d) background. 



Although leptoquark contributions appear in a large number of graphs only a small 
part of them give a sizable contribution to the production cross-section, depending of 
course on the leptoquark parameters such as the mass and the Yukawa coupling. The 
first important consideration is coming from the total width of the leptoquark, which is 
given by 

M Nch 



in 



(12) 



i=i 



where J = 0, 1 is the spin of the particle, / = \ and f\ = \ and N ch = 1 or 2 depending on 
how many channels are available. As is evident for a leptoquark mass of the (9(200 GeV) 
the typical width is of the order of a few hundreds of MeV. This is mainly due to the 
fact that leptoquarks couple to ordinary fermions in a very restricted way, which is not 
generally true for other types of leptoquark couplings, like those appearing in R-parity 
violating MSSM ||. As a consequence of the narrowness of these states the interference 
terms among resonant and non-resonant contributions are rather suppressed. 

The second remark is that, for leptoquark masses M > a/s/2, the single leptoquark 
production mode becomes dominant as far as the first generation leptoquarks are con- 
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cerned. This proceeds mainly via the t-channel graphs, shown in Fig.[| and their con- 
tribution comes mainly from the collinear electrons (positrons). This can be described 



either by the Weizsacker- Williams |26| approximation or by integrating over the momen- 



tum of the final state electron using the leading logarithmic approximation as described in 



reference |24| . We have implemented both approaches, and checked that the results agree 
within less that 10%, which is rather typical for the leading logarithmic (LL) approxi- 
mation we are essentially employing in both cases. As far as the Weizsacker- Williams 
spectrum is concerned, we have used the following form: 

~l + (l-x) 2 , f s (l-x) 2 \ , 2-x 2(x - 1) 



fj/ e [x, s) 2 



In ~^tz 777 + x\n 



x \ m 1 (2 — x) 2 I x x 



(13) 



It should be mentioned that this type of contributions, which are substantially enhanced 
due to the t-channel photon exchange, are only relevant for the first generation lepto- 
quarks and for the semileptonic reactions, Eq.(|lO|)-Eq. (|Tl"l), under consideration. On the 
other hand, first generation leptoquark contributions to Viv^qq and qqq'q' final states are 
relatively suppressed, due to the absence of the t-channel photon exchange enhancement. 

Finally, as the first graph of Fig.0 has a mass singularity in the quark propagator, 
special care has to be paid in the integration in this specific case. As before the LL- 
approximation scheme has been employed. It is worth to mention that comparing with 
the exact result for the total cross-section fll9"l , taking into account the quark masses^, 



and properly integrating over the photon spectrum, both results agree to less than a 
few percent. Moreover in order to have a quantitative estimate of the QCD effects, due 
to the presence of the light quark mass singularity, we compare the perturbative results 
with the resolved photon contribution. In the latter case the total cross section for the 
single-leptoquark production is given by 

a(e + + e~ — > Sq + X) = J dx dzf^/ e {x } s)f q / 7 (z, MIq)ct(xzs) (14) 

which can be written as 

a(e + + e-^S + X) = ^ (±\ f\ ^f l/e ( x , s)f qh (M 2 LQ /(xs), M 2 LQ ) (15) 

where / 9 / 7 is the structure function of the quark inside the photon. In Fig.|| we show the 
cross section as a function of the leptoquark mass for a/s = 192 GeV, for two different 
choices of the structure function parametrizations compared to the perturbative cal- 
culation. We see that, even in the total cross section, one can safely trust the perturbative 
estimate. Moreover this exercise shows us that an error of the order of 10% to 20% is 
to be expected, whereas the differences between the two schemes are concentrated in the 
mass region near the kinematical threshold. 



5 For the light quarks we take m( Uj< 2) = 300MeV. 
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Figure 3: Total cross section for single-leptoquark production, Sq (u-type) as a function 
of the leptoquark mass, for A/e = |. The solid lines represent the resolved photon con- 
tribution (upper: SAS-G ID, lower: DO-G Setl) whereas points refer to the perturbative 
result with light quark mass 300 MeV. 

3 Results & Discussion 

After completing the description of our calculation, we now come to the physics picture 
emerging from it. In the case of NC-type final states, Eq. (|T0"D , the signal has the typical 
structure of a lepton and a jet balancing each other's transverse momentum with some 
hadronic activity in the forward (backward) region. In the irreducible Standard Model 
background, on the other hand, the angular distribution of the positron (electron) is 
peaking in the backward (forward) direction. In order to suppress as much as possible 
the background, without lowering signal's contribution, we have employed the following 
set of cuts: 

m e _ jet > 140 GeV , 5° < max(6 jh 9 j2 ) < 175° , 20° < 6 e < 160° , E e > 5 GeV . (16) 

Among them the most important ones are the cut on the angle of the observed positron 
(electron) and the cut on the electron-jet invariant mass, m e _j et . In Fig.|] we show the 
invariant mass, m e _j et , distribution for leptoquark mass M = 170 GeV and A/e = 0.5, 
with an integrated luminosity L = 500 pb _1 . 

In the case of CC-type final states, Eq . (|Tl"D , the signal, which is much more spectacular, 
consists only of one very energetic jet and a large amount of missing transverse energy. 
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Figure 4: The m e _j et invariant mass distribution in number of events per bin. 



The only cuts employed in this case are 



where px is the transverse energy of the jet. The second cut is employed in order to 
reduce the main irreducible background contribution, coming from single W production, 
in which case two energetic jets, with an invariant mass 0(Mw), would be present in 
the final state at relatively large angles with respect to the beam. Moreover, as the only 
kinematical information available in this channel is the jet momentum, one cannot fully 
reconstruct the mass of the leptoquark state, due to the missing neutrino energyfj. In 
Fig.|| we show the pt spectrum of the signal as well as that of the irreducible Standard 
Model background. 

In summary, we have presented a calculation of leptoquark contributions to two- and 
four-fermion final states at e + e~ collisions, which has been incorporated in the MC event 
generator ERATO. Taking into account the Standard Model irreducible backgrounds we 
have presented a first study of single leptoquark production, in both its NC or CC channel 
at LEP2 energies. 



4 However the leptoquark mass can still be determined by the endpoint of the pt spectrum. 



max(p T x,p T2 ) > 60 GeV , min(^i,% 2 ) < 5° or > 175° , 



(17) 
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Figure 5: The Pt distribution in number of events per bin. 



Appendix: The program ERATO-LQ 



In this appendix we briefly describe the program ERATO-LQ. The skeleton of the pro- 



gram follows closely the skeleton of the mother program ERATO ||22|| . The distributed code 
contains two main directories, f files where all FORTRAN files are stored and run where 
the input files as well as the makefile resides. There are two main FORTRAN files, lqq.f 
and vud.f corresponding to e + e~ — > e~e + qq and e + e~ — > e~v e ud processes respectively. 
The files algor.f, alice,f and common. f contains all necessary routines to make the 
program run 0. The input files look like: 

7,1,1,0 ! process, scalar or vector 

3 ! iterations 

1 !ISR 

128.07 0.2310309 91.1888 2.4974 80.23 2.033 ! input parameters 

1 ! itotal 

-1 . 0.1 !cmax(cmin) cmas 

192 ! energy 



5 Random number generator and the gamma function r(x) are obtained from the Cern Library, 
CERNLIB. 
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50000 !nev 

145. 0.5 !LQ parameters : mass and 1/e 

10000. !scut = (mass-10)**2 

2000 ! IEV 

1 ! iopt :=0 no optimization 

4000 6000 10000 15000 20000 30000 120000 250000 500001 
100. ! cut on invariant mass l+q= mass-10 

0.3 !mass of light quark 

20. 5. 5. 5. 5. 25. ! angles=degrees : clmax cjmax eelmin ejmin c4max cmas 
1 ! naive 



which is very much the same as the one used by ERATO 
the present code, in addition to those already described in reference p2 



The new variables used in 
are: 



I PRO 

This is used to select the corresponding process: IPR0=5 selects e + e" 
IPR0=6 e + e" — >■ e~e + dd. 



e e^uu, 



• ICH0ICE 

This is used to select the specific leptoquark (charge, decay mode) as described in 
Table [I] and Table |. 

• I SCALAR 

ISCALAR=1 if a scalar leptoquark is selected. 

• IVECT0R 

IVECT0R=1 if a vector leptoquark is selected. 

• RLQMAS 

The mass of the leptoquark. 

• CLQ 

The leptoquark coupling, g^gn), as described in Table [I] and Table normalized 
to the electromagnetic one, i.e, g^ <— 1 means gt = e, e = y/4ira em . 

• QMAS 

This is the value of the light quark mass used as regulator of the IR divergences of 
the amplitude. By default is set equal to 300 MeV. 

The final-state kinematics is represented by the three momenta, P4,(e + ), Ps(u) and 
Pio{d). Variable clmax is the minimum angle allowed between the momentum p^ and the 
beam; cjmax the minimum angle (in degrees) allowed between the hardest jet (^3 or p w ) 
and the beam; eelmin is the minimum energy of p^; ejmin is the minimum energy of 
the hardest jet; c4max is the minimum angular separation between the hardest jet and 
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the lepton; cmas is the minimum invariant mass of the two-jet system. If NITER=3 then 
the output provides cross sections for the signal process (ITER=1) and for the Standard 
Model contribution (ITER=2). In addition, the total cross section for on-shell leptoquark 
production is given (ITER=3), where the decay process of the leptoquark has not been 
taken into account and no phase-space cut has been applied. 

Finally a typical output of the lqq. f code, corresponding to the input described above, 
looks as follows: 

TOTAL SIGMA IN NB ERROR 
LQ ERATO 1.749165382052100E-05 3 . 148577170490561E-07 

LQ: on-shell 3 . 939744840638816E-05 3 . 045651348671309E-07 
SM 1.206926939341498E-04 7 . 367659234150230E-05 

[TOTAL, PASS, FAIL ,ICCUT] : 50000 9951 40049 29398 1152 2645 6854 



In addition a programme called phot . f which calculates the total cross section from 
the resolved photon contribution, Eq. (|T5|) , is provided. All the leptoquark coupling (see 
Table [I] and Table |2|) have been incorporated. A link to the standard parton distribution 
functions library, PDFLIB, is assumed. 
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Table 1: Scalar leptoquark charges, couplings and decay modes, as used by ERAT0-LQ 
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Table 2: Vector leptoquark charges, couplings and decay modes, as used by ERATO-LQ 
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